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ABSTRACT 

In this paper we calculate the potential for a prolate spheroidal distribution as in a dark 
matter halo with a radially varying eccentricity. This is obtained by summing up the 
shell-by-shell contributions of isodensity surfaces, which are taken to be concentric and 
with a common polar axis and with an axis ratio that varies with radius. Interestingly, 
the constancy of potential inside a shell is shown to be a good approximation even 
when the isodensity contours are dissimilar spheroids, as long as the radial variation 
in eccentricity is small as seen in realistic systems. We consider three cases where 
the isodensity contours are more prolate at large radii, or are less prolate, or have a 
constant eccentricity. Other relevant physical quantities like the rotation velocity, the 
net orbital and vertical frequency due to the halo and an exponential disc of finite 
thickness embedded in it are obtained. We apply this to the kinematical origin of 
Galactic warp, and show that a prolate shaped halo is not conducive to making long- 
lived warps - contrary to what has been proposed in the literature. The results for 
a prolate mass distribution with a variable axis ratio obtained are general, and can 
be applied to other astrophysical systems such as prolate bars, for a more realistic 
treatment. 

Key words: galaxies: halos — galaxies: kinematics and dynamics — galaxies: struc- 
ture — methods: analytical 



1 INTRODUCTION 

A galactic disc is believed to be embedded in a dark matter 
halo, but the actual shape of the dark matter halo is not yet 
well-understood. Usually, the halo is taken to be spherical 
for simplicity, or oblate. When a deviation from sphericity 
is considered at all, typically an oblate shape is used (e.g., 
Sackett & Sparke 1990). A prolate halo has been used in 
a few cases, e.g. to study warps in galaxies (Sparke 1994, 
Ideta et al. 2000), and to explain the tidal streaming in the 
Galaxy (Helmi 2004). In these the axis ratio is taken to be 
constant, again for simplicity. 

Recently, it was shown that the dark matter halo in the 
Milky Way is prolate with the vertical-to-planar axis ratio 
monotonically increasing to 2.0 at 24 kpc (Banerjee & Jog 
2011). This was obtained by modeling the observed steeply 
flaring atomic hydrogen gas layer in the outer Galaxy, where 
the gas is supported by pressure against the net gravitational 
field of the disc and the halo. It is then natural to ask what 
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would be the dynamical consequences of such a progressively 
prolate halo. This is the motivation for the present paper. 

N-body simulations of galaxies show that the dark mat- 
ter halo is triaxial. Along the Cartesian cuts (X-Y, or X-Z), 
the halos show a range of value from oblate to prolate with 
a preference for a prolate shape (Bailin & Steinmetz 2005, 
Bett et al. 2007, Vera-Ciro et al. 2011). Further the axis 
ratio varies with radius. Indeed, there is no simple physical 
reason why the shape should be constant with radius. Thus 
a radially variable axis ratio for an ellipsoidal mass distri- 
bution is not just of academic interest but of practical use 
for a study of dark matter halos, and in other astrophysical 
systems as well. 

To study the dynamics of a prolate halo, we need to first 
calculate its potential and other physical quantities. Mathe- 
matically this reduces to the problem of inversion of Poisson 
equation for a given density distribution. This is a rich and 
an old subject, which had attracted the attention of mighty 
mathematicians like MacLaurin, Lagrange, Jacobi, Gauss 
to name a few (e.g. Rohlfs 1977), also see Chandraskekhar 
(1969) for details. Yet, interestingly, the calculation for a 
prolate case with varying eccentricity has to our knowledge 
not been studied in the literature so far. 
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Motivated by this, in this paper we study prolate el- 
Mpsoidal meiss distributions with an eccentricity that varies 
with radius. The elHpsoids arc taken to be concentric and 
coEixial or coplanar. We show that the constancy of potential 
within a shell as in homoeoids is a good approximation even 
when the shell is bounded by dissimilar ellipsoids, for small 
variation in the eccentricity as seen in realistic systems. This 
is an important general result from this paper. 

We calculate the potential for a general density distri- 
bution by adding the shell-by-shell contributions from iso- 
dcnsity surfaces from shells internal and external to a point, 
and from this obtain the radial and planar force equations 
and frequencies. Next, these results are applied to the spe- 
cific case of pseudo-isothermal density distribution for the 
dark matter halo, and an exponential disc of finite thick- 
ness, to study the longevity of the Galactic warp for a disc 
in a prolate halo. We show that the prolate halo docs not re- 
sult in a long-lived warp, contrary to the literature where a 
prolate halo (Ideta et al. 2000) or a halo with a radially vari- 
able shape (Binney & Tremaine 1987) have been suggested 
to explain long-lived warps. 

The calculation of the potential and related physical 
quantities are given in Section 2. In Section 3 this is applied 
to toy mass models so as to study the properties of a prolate 
spheroidal system. Section 4 we apply this to the Milky Way 
and study the effect of the prolate halo on the warp and 
show that the halo does not yield long-lived warps. Section 
5 summarizes our conclusions. 

The reader who is mainly interested in the application 
of a prolate halo to the study of warps, may wish to skip 
the technical material involving potential theory in Sections 
2 and 3, and go directly to Section 4. 



2 POTENTIAL OF A PROLATE ELLIPSOID 
WITH VARYING ECCENTRICITY 

The dark matter halo is usually modeled as a spherical or 
oblate distribution of constant ratio, for simplicity (e.g., 
Sackett & Sparke 1990). In a rare case, an oblate mass dis- 
tribution with radially varying axis ratio has been studied 
(Ryden 1990) in the quadrupolar limit, and extended for 
the prolate case by Buotc & Canizares (1996) in the same 
limit. The general case of a prolate spheroidal mass distri- 
bution, in particular with a varying axis ratio, has not been 
studied in the literature so far. The potential for a general 
triaxial case (with constant axis ratios) is stated in Binney & 
Tremaine (1987) without derivation, and has been derived in 
Chandrasekhar (1969), and by de Zeeuw & Pfenniger (1988). 
Prom this one can derive the potential for a prolate spheroid 
of constant axis ratio, as for example was done by Das & Jog 
(1995) for modeling bars in a study of heating due bars. 

In this paper, wc give the analytical calculation of the 
potential and other relevant physical quantities of a body 
for which the isodensity surfaces are prolate in shape with 
radially varying eccentricity. In this way of formulation of 
the problem, the case of a constant axis ratio comes out 
naturally as a simple case of the general calculation, as 
shown later in this section. First the potential due to a 
prolate spheroidal surface is calculated, following the ap- 
proach developed for an oblate spheroidal surface by Binney 



& Tremaine (1987). Then by summing over the contribution 
by the various shells, the net potential is obtained. 

2.1 Potential of a Prolate Spheroidal Surface 

To calculate the potential of a prolate spheroidal surface, it is 
appropriate to use the prolate spheroidal coordinate system 
{u, V, (f>). However, the galactic disc is usually represented in 
cylindrical coordinates. The prolate spheroidal coordinates 
{u, v) and the cylindrical coordinates {R, z) are related by: 



R = A sinhu sinv 
= A coshu cosv 

where A is a constant (see e.g., Arfken 1970). 



(1) 



The Laplacian of the potential can be expressed in prolate 
spheroidal coordinates as: 



du 



d ^ hvh 
du V hu 



d_ 



dv 

hvhu a$ 
d(j> 



(2) 



where the scale factors are: 
hu = A(sinh'^u + sin^v)^^^ 

hif, = A{coshu + sinv) 



(3) 



Here the density is a function of u only, so we expect the 
potential $(w) to be a function of only the coordinate u. For 

potentials of this class, after substituting the scale factors 
in the eq. (2), = reduces to: 



(sinhu = 
du \ du J 



(4) 
(5) 



Hence either 
$ = $0 
which is a constant, or 

(d^/du) = A csch u (6) 
where A is a. constant. On integrating this equation we find, 
' 1 -|- sech u ^ 



$ = -A/2 



/! + sech u\ 

'Of ;~ ) + 'OfV'o 

\ 1 — sech u ) 



(7) 



For large u, sechu — t- 0, and sechu = A/r so a potential of 
the form (eq. 7) varies as: 

$~ {-A/2)(2sechu + log^o) ' 



-A (- - ZosVo) (8) 



where r is the usual spherical radius. 

Hence, if we set logipa — 0, and A = GMj (A), the 
potential given by eq. (8) tends to zero at infinity like the 
gravitational potential of mass M, as expected, hereM is 
the mass of the surface. Thus we are led to consider. 



GM 



whenn<uo 



(9) 



when M > Wo 



This potential is discontinuous on the spheroid w = uo ■ 
Hence it is the gravitational potential of a shell of material 



© 0000 RAS, MNRAS 000, 000-000 



Prolate halos with radially varying eccentricity 3 



on the surface u = uq. This shell has principal semi-major 
and semi-minor axes of lengths Co = Acoshuo and oo = 
A sinhuo respectively. Hence eccentricity of the shell, e = 
sechuo, thus the constant A = coe. We may rewrite eq. (9) 
as: 



$ : 



2co e 



log [j^) when M < Mo, 



(10) 



when u>uo 



On applying Gauss's theorem to the potential across this 
surface, its surface density (following the approach in Binney 
& Trcrnainc 1987, see eq. 2-68) is derived, to be: 



E = 



(11) 



AttG 

where, n„ is a unit vector normal to the surface, and V^" is 
the gradient of the potential: 

1 



Hence, 
E 



A{sinh'^u + sin^vy/^ du 



M 



47ra2g(l — e^cos'^vy/^ 



(12) 



(13) 



2.2 Potential of Prolate Mass Distribution 

The potential due to a prolate spheroidal surface is obtained 
above (eq.[10]). A prolate spheroidal mass distribution can 
be considered as the combination of such surfaces, with the 
condition that the isodensity surfaces should not intersect 
each other. This is so that the function is not unphysical or 
double-valued. Then we can calculate the potential of such 
distributions using eq. (10). 
Define a quantity: 



a (a) 



a dq 
q da 



(14) 



where a is the semi-minor axis. The constraint that the shells 
do not cross is equivalent to stating that the semi-major axis 

of the outer shell must not be shorter than the semi-major 
cixis of the inner shell, this requires that cr (a) > —1. 

Consider a thin shell of matter whose inner surface is 
a prolate spheroid with a semi-minor axis a and axis ratio 
q, and outer surface is a prolate spheroid with a semi-minor 
axis a + 5a and an axis ratio q + dq. The mass of this prolate 
spheroidal shell, when 6a « a, can be shown to be: 



SM = Aira qp{a) 



1 + 



a (a) 



5a 



(15) 



We next calculate E, the surface density of a prolate 
spheroidal shell with dissimilar surfaces following the pro- 
cedure in Binney & Tremaine (1987, see their eq. 2-71) as: 



E = 



p(a)da 
p{a)5a 

772 



a + 



(16) 
(17) 



where R and z are the cylindrical co-ordinates. On using eq. 



(1), and the discussion following eq. (9), these reduce to: 

R — a sinv (18) 
z = aq cosv (19) 



Hence, the surface density of the shell is given as: 



77^ [l + crcx>s\] 



(20) 



a (1 — e^cos'^vy^ 
Using eq. (15) for the mass of the shell, eq. (20) reduces to 



E = 



5M 



A-KO^q (1 



Nl/2 



1 + a cos V 



1 + (a/3) 



(21) 



Note that for similar surfaces, a = Q (eq. [14]), hence the 
surface density of a shell within two similar surfaces ob- 
tained by taking a projection of the mass on the inner sur- 
face of the shell (eqs.[16] and [21]) is identically equal to 
that obtained using the Gauss's theorem approach (eq.[13]). 
Hence the potential inside a shell bounded by similar sur- 
faces can be taken to be constant (Newton's theorem), or 
there is no force at an internal point due to such a shell. 
This is the standard homoeoidal theory of ellipsoids (e.g., 
Chandrasekhar 1969). 

It is not clear a priori that the results from the ho- 
moeoidal potential theory are applicable for a shell bounded 
by dissimilar surfaces considered here. In fact, in this case, 
the surface density obtained by taking the projection (eq. 
[21]) differs from that obtained earlier by the Gauss's the- 
orem approach (eq.[13]) by a dimensionless scale factor S, 
defined as: 



S = 



1 + a cos V 



1 + (a/3) 



(22) 



This scale factor denotes the deviation from the case of 
constant internal potential. We show in Appendix A that 
the values of S are small for the small variation in eccen- 
tricity with radius considered here. Thus, interestingly the 
homoeoidal potential theory though not rigorously valid, is 
still a good approximation. Hence the potential within a 
shell {u < Mo, see eq.[10]) can be assumed to be constant 
even for the case of a shell bounded by dissimilar spheroids. 
In analogy with a homocoid which is a shell bounded by two 
similar spheroids, we define a heteroid to be a shell bounded 
by two dissimilar spheroids. In this paper we study the grav- 
itational potential and force due to a heteroid. 

Hence, for a field point lying inside the spheroid of semi- 
minor axis a (with a prolate co-ordinate Ua), the contribu- 
tion of potential at that point due to this shell (from eqs. 
[9] and [15]) is given as: 



S^int = -2-kG- 
e 



1 + 



cr(a) 



sechua 



sechua 



p (a) 6a (23) 



From eq.(23) onwards, e = sechUa- 

Similarly if the point lies outside the shell, then 



-2^G« 

e 



1 + 



a (a) 



log 



1 + sechu'a 
1 — sechu'a 



p (a) 5a(24) 



where u'^ is the label of the spheroid that is confocal with 
the spheroid with semi minor axis a and lies outside of it. 

The potential of the entire body is the sum of contribu- 
tions from all the homoeoids which make up the body. The 
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contribution from the outer shells at a point {R, z) is: 
^ 5^ir,t = -2-kG 



a>a' 



POO 




/ - 


-' a 





, / 1 + sechua \ , 



da 



(25) 



The contribution from the inner shells to the potential 
at the point {R, z) is: 



-27rG 



Jo e 



1 + 



a{a) 



X log 



1 + sechu'a 
1 — sechu'a 



p{a) 



where a' is given by: 

2 



/2 



a 



R' + 



(26) 



(27) 



The total potential, is given by the summation of eqs. 
(15) and (16) as: 

* = ^ 5*i„t + ^ 5$ex« (28) 
a>a' o<a' 

In the limit of a constant axis ratio, cr = 0, hence 

the above expression simplifies to the result for the prolate 
spheroidal mass distribution of a constant shape. 

2.3 Related Physical Quantities 

Using the above expression for the potential, we next obtain 
other relevant physical quantities such as the force compo- 
nents along R and z, and the angular and vertical frequencies 
at the mid-plane for a particle in this potential. These are 
applied to the particular problem of warps in Section 4. 



2.3.1 Angular Frequency, Q.h 

By taking the gradient of eq.(18), we can calculate the force 
in horizontal and vertical directions. Then force in the radial 
direction is: 



Fr {R, z) = 2nG 



d_ 
dR 



f 

Jo 



1 + 



cr(a) 



1 ~|~ SGchu^^ 



(29) 



The contribution of force from the outer shells is zero since 

the potential is constant inside such shells. To take the gra- 
dient in the cylindrical coordinate system, we rewrite sechu'a 
in terms of R and z, by using the relation, 

R'^sech'^u'' 



(aqe) = 



2 7 2 // 

z sech u„ 



(30) 



1 — sech^u: 

At the galactic niid-plane {z — 0), eq. (30) reduces to: 

«ecW„' = ^j^^JI^_^y/2 (31) 

Substituting eq. (31) in cq. (29), and taking the derivative 
in the radial direction, gives: 

R 



Fr {R) 



47rG r 2 
-^Jo 



1 + 



(7(o) 



p(o) 



[R2 + a2 (g2-l)]i 



7^rfa(32) 



Notice that at the galactic plane, a' = R. 

Then the square of the rotation velocity (Vi? — FrR) is: 

f-R 



cr(a) 



V' {R) = AtxG i aql + 
Jo 

The square of the angular frequency, {fl^ = Fr/R), is: 



2 

a q 







1 



3 



[R2+a2 (g2_i)] 
-Fr 

p{a 



Z:TTTT7^rf«(33) 



[K^+a^ (g2-l)]V 



T7jrf«(34) 



where the subscript h denotes the halo. 
2.3.2 Vertical Frequency, Uh 

The force in the vertical direction is obtained by taking a 
gradient of eq. (18) along z: 



F,{R,z) = 2itG^ 



On simplifying. 



1 + 



d_ 1 + sechu'i ^ 2 ' "'oz'" 
dz 1 — sechu'a 1 — sech'^u'a 



(36) 



On taking the derivative of eq. (30) in the vertical direction, 
we get: 



dsechu'a 



z (^sechu'a — sech^u'a) 



2z^ sech^u'^ - («2 + iJ2 + ^2 (g2 _ 1)) 



(37) 



Substituting eqs. (36) and (37) into eq. (35), and taking the 
limit 2 — >■ gives: 



F4R,z) =47rGz 



Jo 



1 + 



a(o) 



p(a) 



(ii2-|-o2 (g2-l))3/2 
2 



do (38) 



Then the square of the vertical frequency, z/ = F^/z is 
1^1 {R) = 4nG I aq 



1 + 



cr(a) 



p{a) 



{R^+a^ (g2-l)) 



3/2 



da(39) 



For the sake of completeness, in Appendix B a similar cal- 
culation is given for an oblate spheroid with varying eccen- 
tricity. 



3 APPLICATION FOR VARIOUS MODEL 
MASS DISTRIBUTIONS 

In Section 2, we derived the expressions of potential, forces, 
angular frequency and vertical frequency of prolate shaped 
mass distributions. In order to get an idea of how these 
quantities vary with radius, we now apply those equations 
to some simple toy models of mass distributions. We will 
consider that the density distribution varies isothermally as: 



p(a) = 



Po 



1 + ^ 



(40) 



where po and ac are the core density and core radius respec- 
tively (Binney & Tremaine 1987). 

The isodensity surfaces are taken to be aligned prolate 
spheroids whose axis ratio has the value: 
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Figure 1. Isodensity contours in a plane containing the polar axis 
for a prolate spheroidal mass distribution with an axis ratio that 
is (a) constant, (b) increasing, and (c) decreasing with radius. 



qia')=P{l + -r 



(41) 



where a/uc is the dimensionless semi-minor axis and a&/3 
are constants. Within the constraint a{a/ac) > —1 (eq.[14]), 
we use the vahics of a = —0.1, and 0.1 and value of /3 = 1.5. 
Fig. 1 shows the isodensity contours of these mass distri- 
butions. For Q = we get isodensity contours of constant 
axis ratio (g = 1.5) for all radii (Fig. la). For a = 0.1, the 
spheroids become more prolate at larger radii (from q = 1.5 
to 1.9) (Fig. lb), and for a — —0.1 the isodensity surfaces 
become less prolate (from q = 1.5 to 1.2) (Fig. Ic). 

We next obtain the rotation curves for the three cases 
considered. For a density profile with a constant axis ratio, 
the expression for the square of the rotation velocity Vc^ can 
be obtained analytically by integrating eq. (33) as: 



y2 ^ ^-n-Gpoa^ 



1 + e 

1 - e 



c\i/2 /ft 
-I arctan e 







(42) 



At large radii, this is independent of radius, and in the limit 
of i? ^ oo, the constant terminal velocity square is given 
by: 

2nGpoalq 



(g2 _ 1)1/2 

where c and 6 are defined as: 



2/ 2 
= 0-M 



2/ 2 



1) 



(43) 



(44) 



For the case of varying axis ratio, rotation velocity is calcu- 
lated numerically. Fig. 2 shows the rotation velocity curve 
for q: = 0, +0.1 and —0.1. The % differences in resulting 
velocities at large radii for the various values of a are not 
very large. In Section 4, we will apply these results to study 
the pattern speed and winding of the Galactic warp. Here 
we calculate the pattern speed of warp due to dark matter 
halo alone. In the kinematical picture of warps, the pattern 




Figure 2. The rotation velocity for a pseudo-isothermal prolate- 
shaped halo density distribution with an axis ratio that is: con- 
stant (g = 1.5), increasing (q = 1.5 to 1.9) and decreasing (q = 1.5 
to 1.2) with radius. 



speed for the slow mode is given by Q. — u where and v 
represent the angular speed and the vertical frequency re- 
spectively (e.g., Binney & Tremaine 1987). This is obtained 
by treating the warp as a small perturbation and keeping 
only the slow mode since the fast mode will wind up faster. 
The pattern speed of a galactic warp due to dark matter 
halo alone, ((f2p)h), is defined by: 



(45) 



where the subscript h represents the halo contribution. Us- 
ing the density profile and the axis ratio of the dark matter 
halo, as a function of the semi-minor axis from eqs. (40) and 
(41), and substituting these in eqs. (34) and (39), gives the 

square of the angular and vertical frequencies. For the case 
of a constant axis ratio, il.^' = Vc^ / where Vc^ is as given 
by eq. (42), and the vertical frequency square is given by 
integrating eq. (39) to be: 



2 



47rGpoac 
Ve 



(46) 



where 6 and c are defined in eq. (44). 



For density profiles with a varying axis ratio, O/i and Vh 
are calculated numerically. Fig. 3 shows the pattern speed 
curves of different distributions: the pattern speed decreases 
at different rates with radius and approaches zero at large 
radii. For each axis ratio considered, the mass distribution 
near the mid-plane is different which mainly determines the 
rate of fall of the pattern speed curves. The important re- 
sult to note from Fig. 3 is that for a prolate halo the pattern 
speed is > 0, that is, Q.h > Vh- In contrast, for a flattened 
disc or an oblate halo, the reverse is true (e.g, Binney & 
Tremaine 1987). The net value of the pattern speed is de- 
termined by taking account of the contribution from the 
galactic disc also, as will be discussed in Section 4. 



4 APPLICATION TO WARP IN THE GALAXY 

Here we apply the results obtained in Section 2 to the par- 
ticular case of the Milky Way to study the pattern speed 
and winding in the Galactic warp. The aim is to see what 
effect a progressively prolate halo has on the winding time 
of a warp. 
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Figure 3. The pattern speed for the halo-alone case, for a prolate- 
shaped halo density distribution with an axis ratio that is: con- 
stant {q = 1.5), increasing (g = 1.5 to 1.9) and decreasing (g = 1.5 
to 1.2) with radius. 



In the past it has been suggested that a prolate halo 
(Ideta et al. 2000), or an oblate halo with an eccentricity 
that changes with radius (Binncy & Tremaine 1987) can 
help in making warps long-lived. The severity of the winding 
problem depends on the differential precession rate or on 
how rapidly the pattern speed changes with radius. In the 
kinematical picture of warps, the pattern speed for the slow 
mode is given hy Q. — v (Section 3). Wo need to calculate 
the contribution to this both from the disc and the halo. 
In a disc, the pattern speed is negative since v > Q, for a 
gravitating disc, while for a spherical halo the two are equal, 
and for a prolate halo the reverse is true, namely, v < Q,. 
Thus it can be seen that adding a prolate halo may help 
sustain a warp. Qualitatively this was the motivation for 
trying a prolate halo in the present paper. 



4.1 Vertical and angular frequencies due to 
Galactic disc 

For most dynamical studies, the Galactic disc is considered 

to bo infinitosimally thin for simplicity. However, the vertical 
force is then discontinuous at the galactic plane, hence it is 
not possible to approximate small displacements in the ver- 
tical direction perturbatively as a simple harmonic motion. 
Hence we have to use a three dimensional density distri- 
bution. Here we use the simple double exponential density 
distribution for the Galactic disc which is given by: 

p{R,z)=p{Q,Q)exp{-^^ exp{-^^ (47) 

where p(0, 0) is the central density; Ra and Zd are the radial 
and vertical disc scale lengths respectively. 

The corresponding potential at a point (i?,, z) is given 
by Kuijken & Gilmore (1989), also see Sackett & Sparke 
(1990), and Cuddeford (1993): 



$d {R,z) = -^i^Gp{Q,0)Rlzd 



exp 



[l + {kRd 



,2\3/2 



-dk (48) 



where Jo (a;) is the first order Bessol function of order zero. 

The radial force Fii[R,z) is given by —d^d/dR ,and 
then taking the limit as ^; — ^ 0, we get its value at mid-plane 



to be (also see, Kuijken & Gilmore 1989): 
FR{R,0)=47rGp{0,0)Rlzd 

kJi (kR) 



poo 

Jo (1 



(1 + kzd) (1 + {kRdf) 



2\3/2 



dk (49) 



Similarly, for the force along the vertical direction at the 
mid-plane, we obtain Fz{R, z) = —d^a/dz and take its limit 
at z=0, to get: 



F,(i?, 0) = 47rGp (0, 0) R^z 



kJo (kR) 



{1 + kzd) {l + (kRdf) 



2\3/2 



dk (50) 



From these equations of force (eqs. [49] and [50]), the angular 
and vertical frequencies are given by the expressions: 

nd"" = Fr{R,0)/R (51) 

and, 

i^d' = F4R,z^0)/z (52) 

4.2 Disc and halo parameters for the Galaxy 

For the Galaxy, the radial disc scale length Rd — 3.2 kpc 
(Mera et al. 1998) and Zd = 0.35 kpc (Binney & Tremaine 
1987) are used for calculating the force components and the 
frequencies given above. The mid-plane volume density of 
the stellar disc at the solar neighborhood is observed to be 
O.O7M0pc"-^ (Binney & Tremaine 1987). Substituting the 
values of Rd, Zd and p at solar neighbourhood in cq. (47), we 
get the value of the central density, p(0, 0) « 0.99 Mq pc~^. 
Three types of halo distributions are considered: 

(i) A pscudo- isothermal spherical distribution as in cq. (40), 
obtained in the Galactic mass model by Mera et al (1998). 
where po = 0.035 Mq pc~^ and Oc = 5fepc. 

(ii) A prolate spheroidal distribution with an increasing axis 
ratio, that is obtained by modeling the steep flaring of HI in 
the outer Galaxy (Banerjee & Jog 2011); where axis ratio, 
q and the halo density p(a) are given respectively by: 

q = 1 -|- ai (a — 9) -|- a2 (a — 9)^when Qkpc < a < 24kpc, (53) 



and. 



p(a) = 



(54) 



where q = 1 when a < 9kpc, and the constants ai = 0.020 
and a2 = 0.003, and po = 0.035 Mq pc^^ and ttc = 5fepc. 
(iii) A prolate shape with a constant axis ratio, q =1.9 : We 
change the shape of the dark matter halo from a sphere to a 
prolate spheroid by keeping the semi-minor axis fixed such 
that the total mass inside does not change. Thus equating 
the mass inside both spheroids, one can find out the value 
of the core radius Uc for the prolate spheroid. The mass 
enclosed by the prolate spheroid of semi minor-axis a is: 



M(< o) = AnGpoO^qia — ac arctan( — )) 



(55) 



Equating this with the mass inside a sphere (i.e. M(q=l)) 
of radius ttc, gives the value of function of q: 



«=(9).arctan^ = l- 



4 — TT 



ac(l) 



ac(g) 



4(g2 - 1)1/2 2 



Mm 



(56) 
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Figure 4. The pattern speed of warp in tiie Galaxy consisting 
of an exponential disc of finite licigiit and dark matter halo of 
different shapes and distributions. Clearly, the halo contribution 
to the pattern speed is negligible. 



Equating the terminal velocities of the prolate spheroid (eq. 
[43]) and of the sphere gives the value of q{a): 



(57) 



po(g) _ ag(l) (l-g^)V^ 1 
Po(l) al{q) q iZo5(i±f) 

where ac(l) ~ 5kpc and po(l) = 0.035 Mq pc~^ for the 
spherical case {q = 1) (Mera et al. 1998). For q = 1.9, 
this gives Oc = 5.98 kpc and po = O.OIQ Mqpc~^ . For an- 
other typical value, q = 1.5, we get Oc = 5.57 kpc and 
po = 0.021 Mqpc~^. 



4.3 Pattern Speed for Warp 

The pattern speed for a galactic disc embedded in a dark 
matter halo is set by adding in quadrature the contributions 
from both, and is defined to be: 



+ Hd) 



/ 2 



2ni1/2 



(58) 



Using the values of the disc and the three cases of halo 

parameters as discussed above, the resulting pattern speed 
curves arc shown in Fig. 4. The disc determines the main 
behaviour of the pattern speed, hence the latter i.s < 0. Sur- 
prisingly, the presence of dark matter halo makes little dif- 
ference to the pattern speed curve. This is in contradiction 
to the expectation from the literature that halo would affect 
the pattern speed significantly (Section 4.5). 



4.4 Winding Time 

The winding time of a warp between the radii i?i and R2 is 
defined in terms of the pattern speeds at these two radii as: 



i^winding 



7r/[nj,(fli) -f2p(J?2)] 



(59) 



In order to quantify how rapidly warp gets wound up, we 
need to calculate the winding time using cq. (58). If the 
Milky Way were to contain only the galactic disc, then the 
winding time of the warp that extends from 3Rd to 8f?d 
is obtained to be ~ 9.2 x 10^ yr. That is, the warp will 
be wound up within the age of the Galaxy (« 10^° yr). If 
the galactic disc is embedded either in a pseudo-isothermal 
spherical halo or in a prolate spheroidal halo, the winding 



time is almost the same (to within a few percent). This shows 
that the dark matter halo with parameters as chosen here 
does not make any difference in the winding time. Hence, 
if the warp is considered as a bending wave in the galactic 
disc, the gravitational field of the dark matter halo does not 
solve the winding problem of the Galactic warp. 



4.5 Comparison with previous work 

In the past, it has been claimed (Binncy &: Trcmainc 1987, 
Ideta ot al. 2000) that the dark matter halo could increase 
the winding time of kinematic bending wave for a particular 
choice of halo parameters. In these, the galactic disc was 
taken to be infinitcsimally thin, with a surface density that 
falls exponentially with radius: 

-R/Rd 



E(i?) = S„e 



(60) 



where Rd is the disc scale length (the same as in cq. [47]) 
and So is the central surface density. 

For such a disc, at radii much larger than the outer edge 
of the disc, the potential can be approximated as a sum of 
a monopolc and a quadrupole term as shown by Binney & 
Tremaine (1987) to be: 



GMa 



1 + 



mliR" -2z^) 
2r^ 



(61) 



where r is the radius in spherical coordinate system, and 
Md is the total mass of the disc {Ma = 2nEoRa'^). The 
horizontal and vertical frequencies are given by: 



o2 GMd 



1 + 



2i?2 



GMd 

i?3 



1 + 



27Ri 



2R^ 



(62) 



We caution that these arc strictly valid only at large radii 
(>> the size of the disc), whereas in the above references, 
these were incorrectly applied to lower radii within the disc 
~ 10 kpc where the warp begins (see, Ideta et al. 2000, 
Fig. 8; and Binney & Tremaine 1987, Fig. 6-30). Both these 
references use a logarithmic potential for the halo while we 
have derived a rigorous expression for the halo potential. 
Fig. 5 shows the results obtained using this approximation 
(dashed line), and that obtained in this paper (eqs. [51]-[52]) 
for a disc of finite thickness (solid line). Here we used Rd=3.2 
kpc (sec Section 4.1) and Eo = 640 Mopc"^ (Binney & 
Tremaine 1987) for plotting the curves. 

In this approximation, u and Q have comparable values 
- we note that this result is incorrect, since it is well-known 
that within a self-gravitating disc » Q. Check that the 
correct treatment for a finite height disc in the present paper 
does satisfy this criterion (see Fig. 4 which shows fl — v < 0), 
which results in a pattern speed that varies rapidly with 
radius and is dominated by the disc contribution (Fig. 4). 

Figure 6 shows results for the pattern speed of warp 
at the galactic plane for a prolate halo of q = 1.9 and the 
correct treatment for the disc as in this paper (the solid 
line), and it also shows the erroneous results obtained using 
the approximation for the disc (cq.[62]) either by itself, or 
with a spherical halo, and the logarithmic halo as in Ideta 
et al. (2000). The last three cases, in particular the one with 
the prolate halo, show an almost constant pattern speed 
throughout the range of radii of interest where the warp 
is seen. This led Ideta et al (2000) to the conclusion that a 
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Figure 5. and i/^ vs. R for a particle in the Galactic disc using 
the correct treatment as in this paper (solid line), and using the 
approximation of an infinitesimally thin disc and potential valid 
at large radii as in Binney & Tremaine (1987) (dashed line). 
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Figure 6. The pattern speed of warp in the Galaxy consisting 
of a disc of finite thickness and a prolate-halo (q=1.9) from the 
present paper (solid line) , and that using the approximation of an 
infinitesimally thin disc and potential valid only at large radii (as 
in Binney & Tremaine 1987, and Ideta et al. 2000) and different 
halo shapes. 



prolate halo results in a Galactic warp that can be sustained 
for more than the Hubble time. To summarise, they arrived 
at this wrong conclusion due to the use of an approximate 
expression for the potential for an infinitesimally thin disc 
as applied to a low radial range where it is not valid. 

The correct treatment for the disc in the present pa- 
per shows that the disc contribution to the pattern speed 
is dominaiit (solid line. Fig. 6, also see Fig. 4), and that 
the prolate dark matter haJo, indeed halo of any shape, has 
very little effect on the pattern speed of the Galactic warp. 
The normalized differential force due to the heteroid shells 



considered is negligible (< 1%, see Appendix A), hence that 
too cannot affect the above result. In retrospect, this is not 
surprising, since the previous work on the origin of warps by 
non-axisymmetric response showed that the radius of onset 
of warps in a galactic disc depends weakly on whether or 
not an oblate halo is included or on its shape, as shown for 
a disc with an exponential vertical distribution (Saha & Jog 
2006) or a general z distribution (Pranav & ,Jog 2010). 

We caution that we have only considered a simple kine- 
matical model of warp, as an application of the potential 
obtained for prolate ellipsoids of varying eccentricity. The 
latter was the main aim of this paper. A complete treatment 
of the dynamics and origin of warps would be much more 
complex and could involve a responsive and triaxial halo 
(e.g., Binney, Jiang & Dutta 1998, Dubinski & Chakrabarti 
2009) and is beyond the scope of this paper. 



5 CONCLUSIONS 

In this paper, wc have given the analytical calculation for 
the potential of a prolate ellipsoidal mass distribution with 
radially varying eccentricity. This is obtained by summing 
up the sholl-by-shcU contributions of isodensity surfaces. We 
define a heteroid to bo a shell that is bounded by two dis- 
similar surfaces. An important result from this paper is that 
the potential within a heteroid is shown to be approximately 
constant as long as the radial variation and the range of val- 
ues of eccentricity spanned arc small as in realistic systems. 
The normalized differential force across such a shell is non- 
zero but small and docs not affect the dynamical stability 
of such figures on a Hubble time. 

The motivation for this work was to see if the progres- 
sively prolate halo for the Galaxy (as in Bancrjcc & Jog 
2011) could help sustain long-term warps. We apply the 
above calculation to the case of Galactic warp, and on using 
the correct treatment for the disc, find that the disc mainly 
determines the warp pattern speed. Thus, the dark matter 
halo has only a marginal effect on the pattern speed, and it 
does not lead to long-lived warps. 

The results for the potential and force components ob- 
tained here are general and could be applied for a more 
realistic treatment of other astrophysical systems such as 
prolate bars with a varying shape. We also give a similar 
exact calculation for the oblate case for the sate of com- 
pleteness. This can be applied to oblate spheroidal systems 
with varying eccentricity. 

Acknowledgements: 
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important issue about the applicability of the homoeoid po- 
tential theory for the case of spheroids of varying eccentricity 
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Appendix A : Potential and force inside a prolate 
shell within dissimilar surfaces 

In this paper we treat spheroids of radially varying 
eccentricity. In analogy with a homoeoid which is a shell 
bounded by two similar spheroids, we define a heteroid to 
be a shell bounded by two dissimilar spheroids (Section 2.2). 
We show next that the results from the standard homoeoid 
potential theory, namely the constancy of potential inside 
a shell and no net force within it, are good approximations 
even for the case of a shell with dissimilar surfaces - when 
the radial variation in eccentricity is small as seen in realis- 
tic systems. This is an important general result and to our 
knowledge the heteroid case has not been studied systemat- 
ically in the past. 

Potential inside a shell within dissimilar spheroids: Recall 
from Section 2.2 that for the dissimilar prolate surfaces with 
the variation in eccentricity, the dimensionless scale factor S 
(see eq.[21]) denotes the deviation from the case of constant 
internal potential: 



1 + a cos V 

1 + {<^m 



1.06 
1.04 
1.02 
1 



a> 0.98 
re 
u 
to 

0.96 
0.94 
0.92 
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Decreasing q & v=0 

increasing q & v=k/2 

Increasing q & v=0 



123456789 
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Figure 7. The values of the dimensionless scale fa<;tor, S, that 
denote the deviation of the surface density within dissimilaj: sur- 
faces compared to the case of similar surfaces vs. the dimensionless 

major axis, for increasing and decreasing prolatcncss, and for dif- 
ferent values of the prolate co-ordinate v. Note that the deviation 
from the similar shells corresponding to the usual homoeoidal case 
is small < a few %. 



where a denotes the dimensionless change in the eccentricity 
with radius. 

In Fig. Al, we plot the values of this scale factor versus 

the major axis a/oc normalized w.r.t. the core radius, for 
decreasing q (as given by cq. [41]) with a — —0.1 for the 
prolate co-ordinate v= and v = -k (2 and similarly for 
increasing prolateness with q increasing (as given by eq. [41]) 
with a = 0.1. This figure clearly shows that the scale factor 
values are small (< a few %). Hence the deviation from 
the case of similar shells treated in homoeoidal theory is 
non-zero but small. Thus, interestingly, we find that for the 
small variation in eccentricity with radius and for the range 
of q values as considered in this paper- the assumption of 
a constant potential within a shell is a good approximation 
even for the case of a shell bounded by dissimilar surfaces. 

Force mside a shell withm dissimilar spheroids: From the 
homoeoidal theorem wc know that the force inside a shell 
bounded by similar spheroids is zero (e.g., Chandrasekhar 
1969, Binney & Tremaine 1987). Here we find out the force 
due to a shell bounded by two dissimilar spheroids (a het- 
eroid). The expression of the internal potential of homoge- 
nous spheroid is (Binney & Tremaine 1987, eq.2.128): 



#(i^, z) = -nGpila - AiR^ - Asz'^) 



where /, Ai and^3 are defined for the prolate case as follows 
(see Binney & Tremaine 1987, Tables 2-1, 2-2) : 



--'Mm 



1 — e 



_1_ 

2e 



1 + e 



A3 = 2- 



1 + e 



{A3) 



{A5) 



(Al) 



Hence the radial component of the gravitational field gen- 
erated by a homogenous spheroid at an internal point (R,z) 
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Table 1. Normalized differential force magnitudes due to a shell 
within dissimilar surfaces 

Increasing q 



a/ac 


Fr 


Fz 


1 


0.0120 


0.0018 


3 


0.0098 


0.0072 


8 


0.0023 


0.0085 


Decreasing q 






a/ac 






1 


0.0141 


0.0020 


3 


0.0128 


0.0071 


8 


0.0033 


0.0072 



is: 

Fr = -2-KGpAiR (AG) 
The vertical force component is: 

F;. = -2-KGpA3Z {A7) 

We can think of the force at an internal point of a het- 

eroid shell as being obtained by the subtraction of force due 
to the outer spheroid and that due to the inner spheroid. 
Thus, the differential radial force due to the shell is: 

Fr = -2iTGpR{A°i - A\) {A8) 

where the superscripts o and i denote the outer and inner 
spheroids respectively, and where p is the density of the shell. 
Similarly the vertical force is: 

F. = -2-iTGpz{Al - (A9) 

First, check that the force inside of such a shell would 
be identically equal to zero if the surfaces are similar (ho- 
moeoids). Next, we can find out the total differential force 
at a point (R,z) due to outer dissimilar shells by adding the 
contributions of each shell. The magnitude of the differen- 
tial force depends on the eccentricities of ellipsoids which 
surround the shell. The total differential force at a point is 
obtained for the density distribution (eq.[40]) and the eccen- 
tricities (eq.[41]) used in our work, and this is normalised 
with the force due to a homogenous spheroid which has 
the surface that is passing through the point (R,z) and has 
that value of constant density. For normalization, a constant 
density ellipsoid is assumed for simplicity. The increasing q 
(more prolate at larger radii) corresponds to a. = 0.1 while 
decreasing q (less prolate at smaller radii) has a = —0.1 
(see Section 3). These correspond to a variation in the pro- 
late axis ratio, q = 1.5 to 1.9, and 1.5 to 1.2 respectively 
(Section 3). The increasing case covers the range of values 
of q similar to the range deduced from flaring of HI data in 
the outer Milky Way (Banerjcc & Jog 2011, also see Section 
4.2). 

The values for different semi-minor axis are given in 
Table Al. Clearly, the differential force is very small (with 
magnitudes < 1 %) for the range of semi-minor axis values 
that we are interested in. As expected, the differential force 
magnitudes are highest at inner radii but still small ~ 1%. 



The typical value of Oc, the core radius, is 5 kpc in our 
Galaxy (Section 4), thus the ratio a/ac = 1-8 corresponds 
to about 5 to 40 kpc. Thereforc,wc can take the internal force 
due to the spheroids of varying eccentricity to be zero for our 
calculation. This further confirms that Newton's theorem is 
a good assumption even for the case of dissimilar ellipsoids 
considered here. This is an important, general result and can 
be applied to other realistic cases. For example, this is also 
applicable for studies of dark matter halos in simulations 
which show an even smaller variation in the prolate or oblate 
axes ratio (e.g., Jing & Suto 2002, Vera-Ciro et al. 2011). 

Since the typical dynamical timescale in a homogenous 
spheroid with a surface at say 3 flc = 15 kpc is ~ a few x 10® 
yrs, hence the differential force of < 1% of the force due to 
the inner spheroid is not likely to affect the stability of such 
a distribution over 100 times the dynamical timescales or 
over a Hubble time. 

Appendix B : Oblate spehroidal mass distribution 
with radially varying eccentricity 

In this appendix, we give the oblate case with varying 
eccentricity for the sake of completeness, with an approach 
similar to that developed in this paper for the prolate case. 
This problem has been looked at earlier by Ryden (1990), 
who however used the method of multipoles limiting to the 
quadrupole approximations. Thus Ryden's results are ap- 
proximate and in particular not rigorous for inner regions, 
as we check by comparing with the exact calculations below. 
Hence we give the exact calculation below. This could be of 
use not just for dark matter halo but also for other systems 
of varying eccentricity such as bulges. 

The method of derivation of physical quantities for the 
case of oblate spheroidal mass distribution is similar to the 
method that is described for the prolate spheroidal mass 
distributions. So the resulting equations of potential, forces, 
vertical and horizontal frequencies are given here without 
detailed derivations. The density profile, axis ratio and a 
are used as the functions of semi-major axis a. (Remember 
that in the case of prolate distribution, a is the semi-minor 
axis). Here the axis ratio is always less than 1. 

The potential of the oblate spheroidal shell on the sur- 
face u = un in oblate spheroidal coordinate system is given 
in Binney & Tromaino (1987), to be: 

r^m ( arcsin{e) when u < wo, 

$ = _|il£ X I (SI) 

'^''^ [ arcsin (sechu) when m > mq 

where u is defined in cylindrical coordinates as: 

R = A coshu sinv; « = A sinhu cosv iB2) 

where A = aoe is a constant, ao = coshuo, e = sechuo and 
M is the mass on the surface w = wo . 

When the axis ratio varies with radius, the constraint 
that the shells do not cross is again given by cr(a) > — 1 
where cr(a) is defined by eq. (14), except that a is now the 
semi-major axis. 

We next obtain the scale factor for oblate mass distribu- 
tion with varying eccentricity similar to the prolate case (eq. 
[ 22]) to check the validity of the homoeoidal theory. Prom 
above equation (eq.[Bl]), we can calculate the surface den- 
sity of an oblate spheroidal surface using Gauss's theorem. 
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as obtained in Binney & Tremainc 1987, see eq. 2.68): 



(B3) 



Using the equation of spheroid we can find out the expression 

of surface density of a shell bounded by two dissimilar oblate 
spheroidal surfaces by taking the projection of the mass on 
the inner surface of the shell to be: 



p6a 



1/2 



(B4) 



By replacing the cyUndrical co-ordinates by spheroidal 
co-ordinates, we get: 

qp6a 



E = 



^[1 + ' 



(55) 



(1 — e^sin^v) 

Using the equation (15), we rewrite the above equation as 

SMS 



S = 



47ra2 (1 - e^sm^t;)!/^ 



(B6) 



Check that the function S defined here is exactly the 
scale factor as obtained for the prolate case (eq. [22]). On 
comparing eqs. (B3) and (B6), it can be seen that the case 
of dissimilar spheroids deviates from the similar spheroids 
by this scale factor. In this case also, the values of the scale 
factor S would be small. Hence as discussed in the case of 
prolate mass distribution, we can still apply the homocoidal 
potential theory. That is, we can substitute the expression 
of mass of the shell (eq. [15]) in eq. (Bl) to calculate the 
potential due to a shell which is bounded by two dissimilar 
oblate spheroidal surfaces. 

It can then be shown that the contribution of the poten- 
tial at a point (R, z) from outer shells is (this is analogous 
to eq. (26)): 



^ / J a' 

a>a' 



1 H — arcsin{e)p{a)da 

(S7) 

Prom eq. (B7) onwards, e = sechua- 

The contribution from inner shells is (this is analogous 
to eq. (26)): 



a<a' 



r aq 




Jo e 


3 



X arcsin{sechu'^)p(a) da 



where m" is the label of the spheroid that is confocal with the 

spheroid with semi-minor axis a. Wc can replace sech u" in 
terms of the cylindrical coordinates {R, z) using the relation: 



2 l2 // 
, ,2 Z seen Ua , „2 ,2 // 



The potential at {R, z) is 

$ {R, z)=J2 "5^^"* + 



a > a' 

The force in the radial direction is 



a<a' 



<5$e 



(B9) 



(•BIO) 



Fr = AttG 



dR 



f 

Jo 



aq 

e 



l-h 



a(o) 



arcsin{sechu'a)p{a) da 
(SIX) 



At the galactic plane, eq. (B5) reduces to: 



„ ae 
sechUn = — 
E 



(B12) 



and a' becomes R. Substituting eq. (B8) into eq. (B7) and 
taking the derivative gives: 

r-B. 



Fr = 



4tvG_ 
R 



2 

a q 



1 + 



a(a) 



Pi") 



{R-^+a^ {q^-l)y/- 



■ da 



Then the square of the rotation velocity is: 



Vf = 4ivG 



Jo 



1 + 



a{a) 



p{a) 



{R^ + a2 (g2 - l))i/: 



(513) 



■ da 



(514) 

Check that for the special case of a pseudo-isothermal 
density distribution (eq. [40]), and for a constant axis ratio 
so that q{a) = q, the above equation reduces to 



= 47rGpoac9[ 



[i?2 - a2(g2 _ l)jl/2 



(1-92)1/: 



arctan 



Y7^arcsin{l — q 



2x1/2 



- a2(g2 _ l)]l/2 



acQ 



{B15) 

Hence, in the limit R oo, the terminal velocity Vc is: 



V;"(7i^oo) 



AnGpoalq . 2.1/2 /n-,a\ 

(l_g2)i/2 arestn{l -q)' (B16) 



Note that this agrees with the expression obtained by Sack- 
ett & Sparke (1990) for an oblate spheroid of constant shape, 
as expected. 

The square of the angular frequency, = Fr/R is: 



4nG 
52 



2 

a q 



1 + 



cr(a) 



p{a) 



(7?2 + a2 (g2 _ l))l/2 



da 



The force in the vertical direction is 
cr(a) 



(517) 



9zJo ^ 



1 + 



arcsin (^seehu'a) p{a) da 
(518) 

(519) 



Using eq. (B5), we get: 

d sechu'a -lie 
,™ " R{R'^ + a^ (g^-1)) 

Taking the derivative in the vertical direction and using eq. 
(B15), eq. (B14) reduces to: 

F, = 4nGz I 
Jo 



2 

a q 



1 + 



(7(a) 



piu) 



(il2 + o2 (g2 _ 1)) 



3/2 



da 



(520) 



Then the square of the vertical frequency, Vf^ = F^/z is: 



K = 47rG 



I'' 

Jo 



1 + 



a{a) 



p{a) 



{R^ + a^ (g^-l)) 



3/2 



da 



(521) 

The expressions for Fr., Fz, Vc, i^h and are the same as 
those in the case of prolate spheroidal mass distributions. 
The difference lies in the fact that a is semi-major axis and 
the axis ratio q is always < 1. In the prolate case, a is semi- 
minor axis and the axis ratio is always > 1. 
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